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ABSTRACT 

This document provides a user’s guide and program description for the ERTS-B 
Stand-Alone Software Package, This package was developed to compute the 
Keplerian elements at the descending node, the Brouwer mean elements at the 
node, the averaged Brouwer Mean Elements, and other orbital parameters for 
the orbits of the Earth Resources Technology Satellite (ERTS). The typical 
input for the program is an ephemeris file which resulted from a definitive 
orbit generated over a two day arc using the Cowell option in the Definitive 
Orbit Determination System (DODS) or the Goddard Trajectory Determination 
System (GTDS). 
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SECTION 1 - OVERVIEW OF THE ERTS-B 
STAND-ALONE SOFTWARE PACKAGE 


1.1 INTRODUCTION 

The purpose of the Earth Resources Technology Satellite (ERTS-B) Stand-Alone 
Software Support Package is to identify the time of a specific event on a DODS 
(Definitive Orbit Determination System) formatted ephemeris (EPIIEM) file, 
and to compute a set of parameters associated with that event. The specific 
event is defined to be the first descending node which is detected within the 
time span of the file after the epoch. Epoch is defined on the input EPHEM 
file header, or optionally, on cards input by the user. The computed para- 
meters associated with the first descending node include the Brouwer mean 
elements, longitude of the descending node, mean local time, and the 
Brouwer mean elements averaged over a period specified by the user. 

The subroutines used in this package include: 1) DEBTAP (Data Evaluation by 
Trajectory Analysis Program) subroutines (CALSEC/SECCAL, CONVRT, 
CROSSP, DOTPRD, ERRLNG, LAGRIN/ INTRPX, SUBCAP, UNITV); 2) 
modified DEBTAP subroutines (COMPOS, COMPUT, OSCMN); and 3) new 
subroutines (DATOUT, PACK/UNPACK, READER, SHIFT, STORE, YIELD). 
The program was constructed from specifications provided by P. Shapiro and 
E. Herring of GSFC. 

1.2 METHOD 

The ERTS-B Stand-Alone Software Package processing is begun by reading the 
input constants and variables from card. The ephemeris file is then read and 
the first descending node after the specified or default epoch is located. An 
8-point Lagrangiah formula is used to obtain the time of nodal crossing by 
interpolation using the Cartesian z component as the independent variable. 

Next, the Keplerian elements and the related orbital parameters are computed 
at the nodal crossing time. 
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Starting with the point in the EPHEM file closest in time to the descending node, 
a user-specified number of position and velocity vectors, at specified intervals 
of time, are collected. These vectors are converted to Keplerian elements and 
then to Brouwer mean elements and the Brouwer mean semi-major axis, 
inclination, and eccentricity are averaged. A report is then generated of the 
Cartesian, Keplerian, Brouwer mean, and averaged Brouwer mean elements, 
ajong with the related orbital parameters at the time of nodal crossing. A 
functional flow of the program is shown in Figure 1-1, 

1.3 LIMITATIONS 

The ERTS-B Stand-Alone Software Support Package utilizes geocentric true 
of data coordinates. 





Figure 1-1. Program Functional Flow Diagram (1 of 2) 
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SECTION 2 - USER’S GUIDE 


2.1 INTRODUCTION 

The ERTS-B Stand-Alone Software Package requires card and ephemeris file 
input. The input cards are of three types: title card, mandatory constants, 
and optional constants. The title card is used to supply the header on the out- 
put report and may contain any information desired by the user. The mandatory 
constants cards supply the program with values that are not defaulted in the pro- 
gram and must be furnished by the user. Optional constants are those with 
default values defined in the program that may be overridden using the third 
type of card input. These constants are optional but three optional constant cards 
must appear in the input stream, even if blank. 

A spacecraft ephemeris file, in EPHEM format, is the final required input. This 
file may be generated using either the Definitive Orbit Determination System 
(DODS) or the Goddard Trajectory Determination System (GTDS). Program 
output is a printer report of the averaged Brower mean elements and other 
orbital parameters. 

2.1.1 Card Input 

2. 1.1.1 Title Card 

Columns Format Description 

1-80 10A8 Header or title for output report 

2. 1.1. 2 Mandatory Constants Card s 

Input of mandatory constants provides data that are not defaulted in the program. 
These cards must be applied by the user. These values are input in NAMELIST* 
format, 

^ 

Reference: IBM FORTRAN IV Language, #GC28-65l5 
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The NAMELIST name is '&INCARD'. The following parameters are in the NAMELIST: 


Parameter 


Description 


REPTIM 


NREVDY 

NUMREV 

NCNT 


Repeat cycle length, i.e. the total time 
required for a spacecraft (S/C) to repeat 
its ground trace. Ground trace is the 
projection on the earth's surface of the 
spacecraft's trajectory using the center 
of the earth as the center of the projec- 
tion. 

Number of S/ C revolutions per day 

The number of S/C revolutions per repeat 
cycle. 

Number of points to be averaged 


DELMIN Time increment (interval) at which points 

are to be averaged. 

BASLON A pre-assigned longitude of the descending 

node used in calculations to determine 
when a S/C maneuver is required. (Mission 
requirements state that the S/C orbit should 
remain within 18.5 kilometers of this base 
longitude at the specified descending node). 

BASMLT (~) Base local mean time. A pre-assigned 

i - 1, 3 local mean time of the descending node used 

to indicate when a S/C maneuver is 
required, 
i = 1 hour 
i = 2 minute 
i - 3 second 


Units 

days 


rev/daj' 
rev/r. c. 


min. 


deg. 


hr. 

min. 

sec. 


2. 1.1. 3 Optional Constants Cards 


Optional constants are those that have default values defined in the program. 

These constants are input on specific cards and all optional constant cards must 
appear in the prescribed order. A zero or blank field will cause the program to use 
the default value of that constant. There are three optional constant cards, all of 
which must appear in the input deck. 
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Card 1 


Card 2 


Card 3 


Columns 

Format 

Description 

Units 

1-10 

G10 

Epoch (YYMMDD) 

Year, month, day 

11-20 

G10 

Epoch (HHMMSS) 

Hour, min., sec. 

21-30 

G10 

Intermediate output 
indicator 
= 0 no 
* 0 yes 


1-25 

G25.16 

1 

Earth gravity constant 

km 2 /sec 3 

26-50 

G25.16 

1 

Earth radius 

km 

51-75 

G25.16 

2 

Earth dynamic radius 

km 

Columns 

Format 

Description 

Units 

1-25 

G25.16 

2 

Second harmonic constant 


26-50 

G25.16 

2 

Fourth harmonic constant 


51-75 

G25.16 

g 

Earth rotational velocity 

rad. /sec. 


Multiple cases may be processed by stacking the input decks (repeating the above 
sequence of cards for each case in a single execution of the program). Figure 2-1 
shows a sample deck set-up. 


2.1.2 Ephemeris File Input 

The spacecraft ephemeris file input is generated in the EPHEM format by either 
the DODS or the GTDS programs. Figures 2-2 and 2-3 illustrate that format on 
the file. 


■*>- 

HP 


1 

Used for computing Keplerian elements from Cartesian elements. 

2 

Used for computing elements in secular terms which then lead to 
Latus Rectum and its derivatives. 

3 

Used for computing the Greenwich hour angle and the Latus Rectum 
and its derivatives. 
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2.2 OUTPUT 


Program output is a one page printer report of the Cartesian, Keplerian, and 
Brouwer mean orbital elements at the first descending node after epoch and the 
averaged Brouwer mean elements (semi-major axis, eccentricity, and inclination) 
and other mean elements. An example of this report is given in Figure 2-4. 

2.3 JOB CONTROL LANGUAGE 

The card input is read on FORTRAN Reference Number (FRN) 5, the input ephemeris 
file on FRN 24, and the output report on FRN 6. A scratch disk data set must also 
be provided on FRN 30. The program can be executed with the following' JCL: 

// EXEC LOADER 

// SYSLIN DD DSN * GRJJJ. ERTSB. LOAD (AVBME), DISP SHR 
//GO.FT24F001 DD UNIT - 2400-9, DISP = (OLD, KEEP), 

// LABEL - (N,BLP), DCB ~ (RECFM - VS, BLKSIZE - 2808), 

// VOL = SER = XXXXXX 

// GO. FT30F001 DD UNIT = DISK, DISP = (.DELETE), SPACE = ( CYL, (1,1)’ 
//GO.DATA5 DD * 

where N is file number on tape and XXXXXX is tape serial number. 


//ZBWWCIST JOB (GR2111311H, T, G00080, 002002), FFF, MSGLEVEL=(1, 1) 
// EXEC LOADER, REGION=15QK,PARM=EP=MA IN' 

//SYSLIN DD DSN=GTFFM.ERTS.LOAD(AVBME),DlSP=SHR 
//GO. FT24I001 DD DSN=ZBWWCQ01, DISP=(OLD,KEEP), 

// UNIT=(2400-9, , DEFER), VOL=SER=34081, 

// LABEL=(1,BLP),DCB=(RECFM=VS,BLKSIZE»2808) 

/ /GO. FT24F002 DD DSN=ZBWWC002, DISP=(OLD, KEEP), 

// UNIT=(2400-9, , DEFER), VOL=SER =31949, 

// LABEL~(1, BLP), DCB=(RECFM=VS, BLKSIZE-2808) 

/ /GO. FT24F003 DD DSN=ZBWWC003 , DISP=(OLD, KEEP) , 

// UNIT=(2400-9, , DEFER), VOL=SER=35002, 

// LABEL=(l r BLP), DCB=(RECFM=VS, BLKSIZE-2808) 

//GO. FT30F001 DD UNIT=DISK, DISP=1, DELETE), SPA CE=(CYL, (1,1)) 
//GO. DATA 5 DD * 

***** ERTS-B PROGRAM FOR USING EPHEMER1S TAPE 34081 

& IN CARD 


NU MR E V~2 5 1 , REPTIM=18.0, NCNT-20, DELMIN=5, 0, BA SLON=98. 030596, 

BASMLT--9. 0, 42.0,. 15.805 

&END 

750621 200000 0 

3. 986012D+5 6. 378165D+3 

1. 0823D-3 -1.8D-6 

***** ERTS-B PROGRAM FOR USING EPHEMERIS TAPE 31949 

& INCARD 

NUMREV=251, REPTIM=18.0, NCNT^20, DELMIN=5.0, BASLON=98. 030596, 

BASMLT-9.0, 42.0, 15.805 

&END 

750619 200000 0 

3. 986012D+5 6.378165D+3 

1. 0823D-3 -1.8D-6 

***** ERTS-B PROGRAM FOR USING EPHEMERIS TAPE 31949 

& IN CARD 

NUMREV=251, REPTIM=18. 0, NCNT=20, DELMIN=5.0, BASLON=T>8. 030596 

BASMLT=9. 0, 42.0, 15.805 

&END 

750617 200000 0 

3.986012D+5 6.378165D+3 

1. 0823D-3 -1.8D-6 

/* 


Figure 2-1 Sample Stacked Case Set Up 



>R* ** * 


'fc 't* *f* "I 4 




2-5 



2-6 





INTERNAL 

EPHEM OUTPUT 


ENTRY DESCRIPTION (RECORD 1} 

BYTES 

FORMAT 

UNITS 

FORMAT 

UNITS 

(1) 

TYPE OF TAPE ID (THE WORD "EPHEM" 
IN EBCDIC} 

8 

LF 

- 

EBCDIC 

- 

(2) 

SATELLITE NUMBER (7 DIGITS} 

4 

LI 

- 

LF 

- 

(3) 

UTC, A-1 INDICATOR: 
= 1, A.l 
= 2. UTC 

4 

LI 


LF 



(ITEMS 4. 5. AND 6 GIVE START TIME OF 

THE EPHEMERIS) 






(4) 

DATE l 

8 

LF 

DUT 


LF 

YEAR. MONTH. DAY 

(5) 

DAY COUNT OF YEAR < 




' 

LF 

DAYS FROM JAN. 0 

(6) 

SECONDS OF DAY ( 





LF 

SEC 


(ITEMS 7, 3, AND 9 GIVE END TIME OF 
THE EPHEMERIS) 







(7) 

DATE l 

8 

LF 

DUT 


LF 

YEAR, MONTH. DAY 

(8) 

DAY COUNT OF YEAR \ 




)■ 

LF 

DAYS FROM JAN. 0 

(9) 

SECONDS OF DAY ( 



' J 


LF 

SEC 

(10> 

At, INTERVAL BETWEEN EPHEMERIS POINTS 

8 

LF 

DUT 

LF 

SEC 

iiif 

(26) RUM ID TAKEN FROM MIR COMMENT CARDS 

8 

LF 


EBCDIC 
(UP TO 128 
CHARACTERS) 


{27} 

REFERENCE DAY FOR DODS (570918) 

8 

LF 

CALENDAR 

LF 

YEAR. MONTH. DAY 

i28> 

TYPE OF COORDINATE SYSTEM USED 
FOR OUTPUT: 

= 1 , TRUE EQUATOR AND EQUINOX PLUS 
EQUATION OF EQUINOXES 
= 2, MEAN EQUATOR AND MEAN 
EQUINOX 

= 3, TRUE EQUATOR AND TRUE EQUINOX 

4 

LI 


EBCDIC 
AND LI (1ST 4 
BYTES EBCDIC, 
2ND 4 BYTES IN- 
DICATOR VALUE): 
"TRUE" 1 
"MEAN" 2 
"IN-ER" 3 



Figure 2-2 EPHEM File Title Records (1 of 5) 



(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 

(42) 

(43) 


. ■ ' : ; ■. : : ' ; ' .. ■ . . 


INTERNAL 


EPHEM OUTPUT 

ENTRY DESCRIPTION (RECORD 1) 

BYTES 

FORMAT 

UNITS 

FORMAT 

UNITS 

ORBIT THEORY USED: 
= 1, COWELL 
= 2 , BROUWER 

4 

LI 

- 

EBCDIC 

(“COWELL " OR 
"BROUWER “) 

- 

SOLAR RADIATION CONSTANT 

4 

SF 

Gm/OUL- 

DUT 2 

LF 

Gm/Cm-SEC 2 

SOLAR FLUX NUMBER 

4 

SF 

- 

LF 

- 

Cp, COEFFICIENT OF DRAG 

4 

SF 

■ 

- 

LF 

- 

C R , SHAPE PARAMETER 

8 

LF 

- 

LF 

- 

ATMOSPHERE MODEL ID 

8 

LF 

. 

- 

EBCDIC 
("HAR PRIE") 

- 

CROSS SECTIONAL AREA OF SATELLITE 

8 

LF 

DUL 2 

LF 

Cm 2 

MASS OF SATELLITE 

8 

LF 

Gm 

LF 

Gm 

ZONAL AND TESSERAL HARMONICS 
INDICATOR: 

= 0. USED 
= 1. NOT USED 

(USE SAME INDICATOR SETTINGS FOR ITEMS 
38, 39,40,41 AND 42) 

4 

LI 


LF 


BROUWER COMPLEMENTARY PERTURBATION 
INDICATOR 

4 

LI 

— 

LF 

— 

LUNAR GRAVITATION PERTURBATION 
INDICATOR 

4 

■ 

LI 

- 

LF 

’ 

— 

SOLAR RADIATION PERTURBATION 
INDICATOR 

4 

' 

LI 

- 

LF 

' 

- 

SOLAR GRAVITATION PERTURBATION 
INDICATOR 

4 

LI 

— 

LF 

— 

DRAG PERTURBATION INDICATOR 

4 

LI 

- 

LF 

- 

ELEMENTS EPOCH, t Q 

8 

LF 

DUT 

LF 

DUT 


Figure 2-2 EPHEM File Title Records (2 of 5) 
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INTERNAL 

EPHEM OUTPUT 


ENTRY DESCRIPTION (RECORD 1) 

BYTES 

FORMAT 

UNITS 

FORMAT 

UNITS 

(44) 

YEAR \ 

8 

LF 

DUT 

LF 

YEAR 

(45) 

MONTH 1 




LF 

MONTH 

(46) 

DAY V ELEMENTS 




LF 

DAY 

(47) 

(- EPOCH 






HOUR i 




LF 

HOUR 

(48) 

MINUTE 1 




LF 

MINUTE 

(49) 

SECONDS X 1000 / 




LF 

SLC 


(ITEMS 50 THROUGH 55 ARE BROUWER MEAN 
IF ITEM 29 = "BROUWER" CLASSICAL 
KEPLERIAN IF ITEM 29 = "COWELL") 

: 





(50) 

SEMI-MAJOR AXIS AT 

8 

LF 

DUL 

LF 

Km 

(51) 

ECCENTRICITY AT t 

o 

8 

LF 

- 

LF 

- 

(52) 

INCLINATION AT t 

8 

LF 

RAD 

LF 

RAD 

(53) 

ARGUMENT OF PERIGEE AT t 

o 

8 

LF. 

RAD 

LF 

RAD 

(54) 

RIGHT ASCENSION OF THE ASCENDING 
NODE AT l 

o 

8 

LF 

RAD 

LF 

RAD 

(55) 

MEAN ANOMALY AT t 

o 

8 

LF 

RAD 

LF 

RAD 


(ITEMS 56 THROUGH 71 ARE OSCULATING 

VALUES EVALUATED AT t ) 
o 






(56) 

TRUE ANOMALY 

8 

LF 

RAD 

LF 

RAD 

(57) 

ARGUMENT OF LATITUDE = (53) + (56) 

8 

LF 

RAD 

LF 

RAD 

(58) 

FLIGHT PATH ANGLE 

8 

LF 

RAD 

LF 

RAD 

(59) 

ECCENTRIC ANOMALY 

8 

LF 

RAD 

LF 

RAD 

(60) 

PERIOD 

8 

LF 

DUT 

LF 

DUT 

(61) 

PERIGEE HEIGHT 

8 

LF 

DUL 

LF 

Km 


Figure 2-2 EPHEM File Title Records (3 of 5) 




EPHEM OUTPUT 


UNITS 

FORMAT 

UNITS 

DUL 

LF 

Km 

RAD/DUT 

LF 

RAD/DUT 

RAD/DUT 

LF 

RAD/DUT 

RAD/DUT 

LF 

RAD/DUT 

DUL 

LF 

. . 

DUL 

DUL/DUT 

LF 

DUL/DUT 

DUT 

LF 

DUT 

2 

RAD/DUT 

LF 

RAD/DUT 2 

OR UNITS 


OR INTERNAL 

0F/> ; 


UNITS GF p. 

RAD/DUT 3 

LF 

RAD/DUT 3 

- 

LF 

- 


SIXTY 8-BYTE 
EBCDIC WORDS, 
ONE FOR EACH 
UNKNOWN 
USED) 


DUT 

LF 

DUT 

DUT 

LF 

DUT 



Records (4 of 5) 


I 
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ENTRY DESCRIPTION (RECORD 1) 

INTERNAL 

EPHEM OUTPUT 

BYTES 

FORMAT 

UNITS 

FORMAT 

UNITS 

(195) At, INTERVAL BETWEEN EPHEMERIS POINTS 

8 

LF 

DUT 

LF 

DUT 

(196) LPN INDICATOR 

4 

LI 

— 

LF 

— 

= 1, INCLUDE PRECESSION AND NUTATION 






= 2. INCLUDE PRECESSION ONLY 






= 3, NEGLECT PRECESSION AND NUTATION 






(197) GREENWICH HOUR ANGLE AT EPOCH 

8 

LF 

RAD 

LF 

RAD 

(198)-(293) NOT USED (ZEROS) 

8 

LF 

- 

LF 

- 

/ COSINE IDENTIFICATION 

8 

LF 

- 

) 

- 

1 COSINE HARMONIC COEFFICIENT 

4 

SF 



1 

_ 

(29«)-(296) { 




/ lf 


I SINE IDENTIFICATION 

8 

LF 

- 


- 

' SINE HARMONIC COEFFICIENT 

4 

SF 

— • 

J 

“ . 

/ COSINE IDENTIFICATION 

8 

• 

LF 

• 


- 

\ COSINE HARMONIC COEFFICIENT 

4 

SF 

__ 

| 

_ 

(348M350) ' 




> lf 


I SINE IDENTIFICATION 

8 

LF 

- 


- 

* SINE HARMONIC COEFFICIENT 

4 

SF 

• 

. — • 

) 

— 


ENTRY DESCRIPTION (RECORD 2) 


THIS RECORD IS A CONTINUATION OF RECORD NUMBER 1. THERE IS A MAXIMUM OF 135 PAIRS OF SINES AND COSINES OF 
HARMONIC COEFFICIENTS, OF WHICH THE FIRST 19 APPEAR IN RECORD NUMBER 1. THE REST, THAT IS, 20-135, ARE IN RECORD 2. 
THE FORMAT IS AS SHOWN IN RECORD NUMBER 1. 
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DOUBLE 




INTERNAL 


INPUT/OUTPUT 

WORD 















NUMHI-R 

ENTRY DESCRIPTION (RECORD 3) 

BYTES 

FORMAT 

UNITS 

FORMAT 

UNITS 

1 

DATE OF FIRST EPHEMERIS POINT 

TIME t OF 

js 

LF 

OUT ( 

LF 

YEAR, MONTH. 
DAY 

2 

DAY COUNT OF YEAR j 

FIRST EPHEMERIS 

I 


j 

LF 

DAYS FROM JAN 0 

3 

SECONDS OF DAY 

1 POINT 

\ 



LF 

SEC - 

4 

INTERVAL BETWEEN DATA POINTS (At) 

8 

LF 

DUT 

LF 

SEC 

5 

X ) 


8 

LF 

DUL 

I.F 

DUL 

6 

v ( POSITION VECTOR COMPONENTS AT 
Y 1 TIMEt 

8 

LF 

DUL 

LF 

DUL 

7 

2 ) 


8 

LF 

DUL 

LF 

DUL 

8 

X ) 


8 

LF 

DUL/DUT 

LF 

DUL/DUT 

9 

• ( VELOCITY VECTOR COMPONENTS AT 
: ( TIMEt 

8 

LF 

DUL/DUT 

Lr 

DUL/DUT 

10 

z J 


8 

LF 

DUL/DUT 

LF 

DUL/DUT 

11-304 

49 SETS OF POSITION AND VELOCITY VECTORS AT 

48 PER 

LF 

DUL. 

LF 

DUL, 


TIMES t + At, t + 2At, . . . , t + 49 At 


SET 


DUL/DUT 


DUL/DUT 

305 

TIME t OF FIRST EPHEMERtS POINT ! 

N DUT 

3 

• 

LF 

DUT 

LF 

DUT 

306 

I At |. TIME DIFFERENCE BETWEEN POINTS IN DUT 

0 

LF 

DUT 

LF 

DUT 

307-350 

ZEROS (NOT USED) 


3 

— 

- 

— 

- 




EACH 







Figure 2-3 EPHEM File Data Record 








ERTS-B STANO-ALONE SOFTKAWd PACKAGE FOR USING EPHEMEBIS TAPE 35002 - — RECORD 1 
SATELLITE ID. =7205301 



TAPE 

BEGIN 

X 

JUN. 

16. 

1975 

0 

HR. 

0 

MIN. 

0.000000 

SEC. 


TAPE 

END 

= 

JUN. 

18 . 

1975 

6 

HR. 

0 

MIN . 

0.000000 

SEC . 


tape 

EPOCH 


JUN. 

17. 

1 975 

20 

HR. 

0 

MIN . 

0.000000 

SEC. 

INPUT 

SEARCH 

TIME 

• — 

JUN. 

17. 

1975 

20 

HR. 

0 

MIN . 

0.000000 

SEC. 

AFTER 

SEARCH 

TIME 

•3£ 

JUN. 

17 . 

1975 

20 

HR. 

58 

MIN. 

31.292060 

SEC . 


GEOCENTRIC TRUE OF DATE COORDINATES 


to 

i 

M 

to 


K= 0.5 1017376425621510 
XDQT= 0.32422294381 204670 


TYPE OF ORBIT = ELLIPTICAL ORBIT 


0* Y= 0.52044420 10535604D 04 

00 YOOT- -0.005678598 1227023D 00 


Z» 0.0 

ZDOT= — 0 .7308473334 05722 9D 01 


HT= 0.909791940 OJ LAT = 

SMA= 0.729446410 04 =C C= 

RCA= 0.7 Za 7 7 3 760 04 APR= 

HCA= 0.909597760 OJ APH = 


0.0 

0.922I225ID-03 
0. 7 JO 1 1 9050 04 
C.92J050S40 03 


LDN= 0. 18 54433 3D 03 
AJ NC= 0. 96 96205 90 02 
PER = 0. 17 22230 ’D 01 
HCAN= 0.49114J50D 03 


PREDI CTFD 
PREDICTED 


T l ME 
TIME 


OF 

UF 


V = 0 .739800020 01 GAM = 

LAN = 0.225570960 03 APF= 

TA= 0.138 132 1SD 02 EA= 

APHN = 0 *498 407 4 2D 03 


0.1 2603 1600-01 AZ~ 
0.1 66 I 8678D 03 SML= 
0.13800607D 02 M4= 


PER I CENTER PASSAGE 
PE R I CENTER PASSAGE 


BEFORE DESCENDING NODE 
AFTER DESCENDING NODE 


17: JUN: 75 
17:jun:7s 


0. 1 38962050 
0. 729445790 
0«- *37880030 


20 

22 


54 

37 


9.330270 

0.426300 


SMA = 0.728558140 04 KM 

ECC = 0.617641150-03 

INC = 0.969677533 02 UtG 

L. AN =• 0.225570950 03 Ucu 

APF = 0.156198250 03 OcU 

MA = 0.217765850 02 OlG 

APGH = 0.911941260 03 KM 

ERRLON 


— 0.4 JJ0 10230 01 N MILES 


NOOE RATE 

PERIGEE PATE 
ANOMALY RATE 
MEAN MOTION 
OVERLAP 
RPTCYL ERROR 
PRGH 
E RRMLT 


ME AN 

0.97231 305D 
'0.274 106070 
0.50229 4770 
0.101 524940 
0.150005000 
0. 1 3.J 15 4750 
1.902R4 1510 
-0.531 107 75D 


.LMENTS 

00 OEGVDAY 

01 DEC./-OAY 
04 OE 020 AY 

-02 WAD/SEC 

02 PERCENT 
01 N MILES 

03 KM 

03 sec 


DLR = -0.258166310 02 DEGZREV 
DLROA = -0.534077620-02 DEG/’REV/KH 
CLRDI * 0.84035691D-02 DEG/’REV/'DEG 

DOMGDA » -0.334429550-04 DEG/’REVZKM 
OUMGD1 = 0. 771520410-02 DEG/REV/OE6 

MEAN LOCAL TIME •= 9 20 O.t0892O 02 

NODAL PERIOD ^ 0.I0326270D 03 MIN 

ANOMALISTIC PERIOD * 0.10314632D 03 MIN 


Figure 2-4 Sample Program Output (1 of 3) 
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ERTS-d 5 T AND— ALDNt SOFTWARE PACKAGE FOR USING EPHEME R IS TAPE 31949 


RECORD 


SATELLITE ID. =* 7205e01 



TAPE 

dEGlN 

X 

JUN. 

IQ. 

1975 

0 

HR. 

0 

M IN. 

0.000000 

SEC 


TAPE 

END 

X 

JUN. 

20 . 

1 975 

6 

HR. 

0 

MIN. 

0.000000 

SEC 


TAPE 

EPOCH 

X 

JUN. 

19 . 

1975 

20 

HR. 

0 

MIN. 

0.000000 

SEC 

I NPUT 

SEARCH 

Tl ME 

= 

JUN. 

19 . 

1975 

20 

HR. 

0 

MIN. 

0.000000 

SEC 

AFTER 

SEARCH 

TIME 

= 

JUN. 

19 . 

1975 

2 1 

HR. 

9 

MIN. 

52.632519 

SEC 


GEOCENTRIC TRUE OF DATE COORDINATES 


M 


l-i 

CO 


TYPE QF OR3I T 


ELLIPTICAL 0R8IT 


X= 0.4921022736330473D 04 
XOOT= O .a5151Y7di9taj976D 00 


Y — 0.5374770771 22 70 2 50 04 

Y JOT = -0.77692400870001390 00 


Z~ 

2DOT* 


0.0 

-0.730911801 85070850 


01 


HT* 0.9091 52 0 63 03 
SHA= 0.729447140 04 
RCAF 0.7 28 7 02 b6 J 04 
HCA = 0.908886790 03 


LAT= 0.0 

LCC= 0 . 1 0 20 58 1 2D' 
APR— C- 7301 31 60D 
APH= 0.923775990 


LON = 0. 1 32 58204 D 03 V= 0.73994529D 01 
■02 AINC = 0.989620573 02 LAN= 0.227523460 03 
04 PER= 0.17 222 6 1 2D 01 T A= 0.15356E77D 02 
03 HCAN= 0.49075960D 03 APFN = 0.4987S913D 03 

PRcDICTfcO TIME UF PER [CENTER PASSAGE BEFORE DESCENDING 
PREDICTED TIME OF PER [ CENTER PASSAGE AFTER DESCENDING 


GAM= 0. 1 54707330-0! 
APF= 0. 1646431 2D 03 
EA= 0.153413990 02 


AZ= 0.18896206D 
SML— ’ 0.729446330 
MA= 0.1 53259280 


NODE •= 
NODE « 


19: JUN 175 21 5 0.28S80D 

19:jUN:75 22 48 0.37694D 


03 

04 

02 

02 

02 


MEAN ELEMENTS 


SMA * 

ECC = 
INC - 
L AN - 
AP= « 
MA = 
APGH - 
: RRLUN = 


0.720553470 04 

0.649327613-33 
0.989670673 02 
0. 22752 3453 03 
0. 171377710 0 3 
0. 059 13 3293 01 
0. 91 2 l 754 ID 03 
-0. 4203964 1 D 01 


KM 

DEG 
OcG 
JE a 
Udu 
KM 

N -MILES 


NODE RATE 
PERIGEE RATE 
ANOMALY RATE 
MEAN MOTION 
OVERLAP 
WPTCVL ERROR 
PRGH 
l RR MLT 


0.37 230 23 H D 00 
-3.274197280 01 
.0.502294430 04 
3. 101 52 4 88 0— 0 2 
3.150053240 02 
0. I. 2559 1 1 20 0 1 
0.9 027 1 394 D 03 
-0. 132328200 04 


DEG/ 0 AY 

DEG/DAY 

DEG/DAY 

RAD/SEC 

PERCENT 

N MILES 

KM 

SEC 


-0.250166500 
— 0.534 C777 2D 
0.84 03558 OD 
-0.334425710 


DLR 
OLRDA 
, CLROI 

COM G DA 
DOMGDI 
MEAN LOCAL TIME = 
NODAL PERIOD 
ANOMALISTIC PERIOD =* 


= 0.771519420-02 


OEG/RE V 
OEG/REVKKH 
DEG /REV/ DEG 
DEG/REV/KM 
0EG/REV2DEG 
9 20 0.12523D 02 
0.103262770 03 MIN 
0.103146890 03 Ml N 


02 

02 

•02 

04 


Figure 2-4 Sample Program Output (2 of 3) 
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ERT5-0 it AND-ALONE SOFTWARE PACKAGE FOR USING EPHEMESIS TAPE 34081 — — RECORD 1 


SATELLITE ID. = 7205BO 1 


TAPE BEGIN = JUN. 20. 1975 0 hR. 0 MIN. 0.000000 SEC. 

TAPE END - JUN. 22. 1975 6 NR. 0 MIN. 0.000000 SEC. 

TAPE EPOCH - JUN. 21. 1975 20 NR. 0 MIN. 0.000000 SEC. 

INPUT SEARCH TIME = JUN. 21. 1975 20 NR. 0 MIN. 0,000000 SEC. 

FIRST DESCENDING NODE AFTER SEARCH TIME = JUN. 21, 1975 21 HR. 21 MIN. 14.445623 SEC. 


GEOCENTRIC TRUE OF OATE COORDINATES 


♦ * * * 


TYPE OF ORUIT s ELLIPTICAL ORBIT 


X= 0.47 346341741 60552D 04 Y = 0 . 55 309 144 I 032 1 4 4 4D 04 Z= 0.0 

XDQT = 0.Q77P4JI 4741 5SJ69D 00 YOOT= - 0 . 74 7 1 506300 0 8 573 2D 00 ZDOT= - 0.T3096955I 6236487D 01 


HT= 0. 906583100 03 LAT= 0.0 • LGN= 0.1797159ID 03 V= 0.74000332D 01 GAM*= 0. 1 898 04 l ID- 0 1 AZ= 0. 18B961Q4D 0} 
S MA — 0.72944776D 04 ECC= 0 . 1 1 l 34 7 92D- 02 AINC= O.9R961d30D 0? LAN= 0.22947636D 03 APF= 0.I6267294D 03 SML= 0.72944685D 04 
RCA= 0.7 26 635530 04 APR= C.7J025998D 04 PER = 0.I7222634D 01 0.17327060D 02 EA= 0.173030700 02 -MA^ 0.172390390-02 
HCA= 0 .9032 1 5 J2u 03 APH= 0.924459B2D 03 HCAN= 0.4RJ39T04D 03 APHN= 0.499 16837D 03 

PRE01CTE0 TIME OF PER I CENTER PASSAGE BEFORE DESCENDING NODE = 2i:JUN:75 21 16 0. 16602D 02 
PRtOICT-D TIME OF PERI CENTER PASSAGE AFTER DESCENDING NODE - 21:JUN:7S 22 39 0. 234960 02 


* 4 4 4 4 4 4 4 4 * « A 4 4 4 444 4 4 44 4444*44444 4 444 4*444444 


MEAN elements 

SMA = 0.72855654D 04 KM NODE RATE = 0.9 72260760 00 DEG/DAY DLR = -0.258166560 02 DEG/REV 

ECC = 0 .&•, 34 1 4030-03 PER IGEE RATE = -0.274 19961D 01 DEG/CAY ' CLRDA = -0 .53407768D- 02 DEG/REV/KM 

INC - 0 • 5d V 6735 70 02 Duo ANOMALY RATE = 0.5J2294J5D 04 DEG/DAY CLRDI = 0.3 4 035*0 30-02 DEG/REV/DEG 

LAN = 0.22947&35J 0 3 Dc-u. MEAN MOTION = 0. I 01 52 4860-02 HAD/SEC DOMGDA = - 0 . 334 4 1 4 1 30- 0 4 DEG/REv/KM 

AP= = 0. 170526350 03 Dc-. OVERLAP = 0.15059924D 02 PERCENT DOMCDI = 0 . 77 1 S 1 9 870- 0 2 DEG/REV/OEG 

MA = 0.943524320 01 0£C RPTCYL ERROR — 9.1I5921S9D 01 N MILES MEAN LOCAL TIME = 9 20 0.62400D 01 

APGH - 0.912497323 03 KM PRGH = 3.302393460 03 KM NODAL PERIOD = 0.103262790 03 MIN 

ERRLUN = -0. H0662372D 01 N MILES ERRMLT = — 0 » l 3295«i420 04 SEC ANOMALISTIC PERIOD = 0.103146910 03 MIN 


Figure 2-4 Sample Program Output (3 of ;>} 


SECTION 3 - PROGRAM DESCRIPTION 


The ERTS-B Stand-Alone Software Package consists of an executive routine, 
fifteen subroutines, one function subprogram, and one BLOCK DATA subprogram. 
The Executive, the BLOCK DATA, and subroutines DATOUT, YIELD, PACK, 
READER, SHIFT, and STORE were designed and programmed by the Computer 
Sciences Corporation (CSC); subroutines COMPUT, COMPOS, and OSCMN 
were extracted from the DEBTAP subroutines of the same names; and sub- 
routines CALSEC, CONVRT, CROSSP, ERRLNG, LAGRIN, SUBCAP, UNITV, 
and function DOTPRD were taken from DEBTAP with no modifications. Only 
the CSC subroutines and the modified DEBTAP subroutines are documented 
herein. The unmodified DEBTAP subroutine descriptions may be obtained 
from the Goddard Space Flight Center, mail code 571. 2. 

3.1 SUBROUTINE DESCRIPTIONS 

The subroutines that were designed and programmed or modified by CSC are 
described in the following pages . 
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Subroutine: Executive Routine 

PURPOSE 

« . 

To compute the averaged Brouwer mean elements at the first descending node. 

METHOD 

The main routine is divided into four parts. Each part has a specific function. 

Part I: a) Reads in data from both card deck and ephemeris tape. 

b) Searches for descending node from epoch or tape start whichever 
occurs last. 

c) Uses Lagrangian interpolation to obtain nodal elements by assuming that 
the Z coordinate is the independent variable and setting Z equal, to zero. 

d) Obtains the closest point to the nodal point according to time. This 
then becomes the first point used for averaging. A number of sub- 
sequent points are then collected over a given interval. 

Part II: a) Calls subroutine COMPUT to obtain Keplerian elements from Car- 

tesian elements at the descending node, 
b) Calls subroutine OSCMN to obtain averaged Brouwer mean elements 
from those points collected previously. 

Part III: Collects and outputs appropriate parameters on a single page 

Part IV: Processes error conditions. Should any error occur an appropriate message is 

printed and the run terminated. Should no error occur the remainder 
of the ephemeris file is read and discarded and the control is then trans- 
ferred to the beginning of the program to start another ephemeris file. 
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USAGE 


X. Calling sequence: 

Not applicable 

2. COMMON blocks used: 

$BLANKCOM, CBSTAT, CNVFCT, DUMP, ECONS, ELEM, 
INDATA, HAS 

3. COMMON block parameters used: 


FORTRAN 

Name 

COMMON 

Name 

I/O 

Dimension Description 

KTR 

CBSTAT 

I/O 

400 

Subroutine/entry point counter 

DEGPRD 

CNVFCT 

I 

1 

Degrees per radian 

TWOPI 

CONFCT 

I 

1 

Pi times two 

L 

DUMP 

0 

1 

Intermediate output control 
30 = no intermediate output 
6 = intermediate output 

IRAD 

ECONS 

I 

1 

Radius of Earth in km 

EMU 

ECONS 

I 

1 

Gravity constant of Earth in 
km^/sec 

X 

ELEM 

I/O 

100 

I Position vector of a point in 

j km 

Y 

ELEM 

I/O 

100 

Z 

ELEM 

I/O 

100 

XDOT 

ELEM 

I/O 

100) 

Velocity vector of a point in 
km/sec. 

YDOT 

ELEM 

I/O 

100 

ZDOT 

ELEM 

I/O 

100 1 

TITLE 

INDATA 

0 

10 

Title record 

DELMIN 

INDATA 

0 

1 

Tentative interval size between 
point in minutes 

NCNT 

INDATA 

0 

1 

Number of points to be used for 
averaging 

IEPYMD 

INDATA 

0 

1] 

Optional epoch in YYMMDD and 

IEPHMS 

indata 

0 

ll 

HHMMSS 

I WRITE 

INDATA 

0 

1 

3-3 

Intermediate output indicator that 
defines L in common block DUMP 
I WRITE = 0, L - 30 
no intermediate output 
I WRITE * 0, L = 6 
intermediate output 



4. Subroutines and functions used: 


CALSE C/SEC CAL, COMPOS, CROSSP, DATOUT, 

DOTPRD, ERRLNG, LAGRIN/lNTRPX, OSCMN, PACK/ 
UNPACK, READER, SHIFT, STORE, SUBCAP, UNITV 

5. This routine is called by: 

Not applicable 

6. Core locations required: 

11984 bytes 

COMMENTS AND RESTRICTIONS 

The output statements are partially taken from subroutine CLASEM of the DEBTAP 
program. 

REFERENCES 

Documentation and program listing of subroutine CLASEM of the DEBTAP program. 












Subroutine: Compute Keplerian Elements and Other Related Parameters 
(COMPOS) 


PURPOSE 

To compute Keplerian elements and other related parameters from Cartesian 
elements. 


METHOD 

The Keplerian elements and other related parameters are computed as follows: 
1. Semimajor axis (SMA) 

position 


r = Jr • r 
o v o o 


/ • • 
= v r • r 


o o 


r • r 

D = -2 — 

0 yiT 


velocity 

parabolic anomaly 


o 

U 


r . r 
o o 

U 


a = 


2 - r 


o 

U 


= SMA 


where U is the gravitational coefficient of the Earth. 
2. Eccentricity (ESQRT) 


e 2 = (l-J£ 


e = ESQRT 


+ — = ECC = ESQRT 2 
a 


3. Semiparameter (SEMIP) 

P=r o( 2 -T)- D o =Ml-e 2 ) 

4. Point of closest approach (RCA) 



I D 2 

2 o 


P 

1 + e 


= RCA 


HCA = 
HCAN = 


q - RADIUS 

HCA 

1.852 


5. Apocenter radius (APR) 


q 0 = a (1 + e) = APR 

« i 

APH = q - RADIUS 


APHN = 


APH 

1.852’ 


6. True anomaly (TA) 



1 sin f 

f = tan" 1 — * TA 

0 cos f 

o 


= SEMIP 
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7. Orthogonal vectors (P and Q) 

The P and Q vectors are formed from the orthogonal unit vectors U and V 
and the true anomaly by: * 


P = U cos f - V sin f 
o o 


Q = U sin f + V cos f 
o o 


where 




UI+UJ+UK 
x y z 


and 


V 


rr - rr 




= (jZ02l)t + (sLliAj + (ihr J*. 
\ Vpu 7 \ Vpu / \ Vpu 


K 


8. Inclination (AINC) 
-1 


i = tan 


v 


U 2 + V 2 
z z 


7 


(U + V ) 2 + (U - V ) 2 >1 

x y y x 


= AINC 


9. Longitude of the ascending node (ELAN) 

The true longitude is defined by: 

JH =£2 + co + f = S2 + u 
o o o 

where 

fi = longitude of the ascending node 
oj = argument of perigee 
u = argument of latitude 
f = true anomaly (computed in step 6 as TA) 


3-8 



True longitude is computed in terms of the orthogonal vectors U and V as: 


Jl - tan 
o 


- / U - V 

-i f y x 

V u +v 

NX y 


The argument of latitude is: 


u = a) + f = tan 
o o 




thus yielding 


52 = 1 - u = ELAN 
o o 


10, Argument of perigee (APF) 


co = u - f = APF 
o o 


11. Flight path angle (GAM) 

/ . ' . v , / e sin f 

Y = tan" 1 ( ■ — ~ ) = tan" 1 ' 0 


cos Y 


1 + e cos f 


= GAM 


12. Azimuth (AZ) 
-1 


A = tan 


r Q (xy - xy) 


[y (zy - zy) - x (xz - xz) 


- tan 


-1 


r (rW ) 
o z 


y(rW x )-x(rW y ) 


where W , W , and W are the components of W, the cross product of r and r 
x y z o o 

The quadrant of the azimuth is checked by the library function DATAN2 as follows: 


n (4) 


(t) 


m 




IV 


(t) 


13. Period (PER) 


P = 




2 a 3 


= PER 


3600 



USAGE 


1. Calling sequence: 

CALL COMPOS (RVEC, VVEC, ELM) 


FORTRAN 

Name 

I/O 

Dimension 


Description 

RVEC 

I 

3 


Position vector 

VVEC 

I 

3 


Velocity vector 

ELM 

o 

6 


Keplerian elements 

2. COMMON blocks used: 



CBSTAT, CNVFCT, DUMP, 

ICONS 


3. COMMON block parameters used: 


FORTRAN 

Name 

COMMON 

Name 

I/O Dimension 

Description 

KTR 

CBSTAT 

I/O 

400 

Subroutine/entry point counter 

DEGPRD 

CNVFCT 

I 

1 

Degrees per radian 

TWOPI 

CNVFCT 

I 

1 

Pi times two 

L 

DUMP 

I 

1 

Intermediate output control 
30 = no intermediate output 
6 - intermediate output 

ERAD 

ECON 

I 

1 

Radius of Earth in km 

IMU 

ECON 

I 

1 

Gravity constant of Earth in 
km^/sec^ 


4. Subroutine and functions used: 

CROSSP, DOT PRD, UNITV 

5. This routine is called by: 


OSCMN 

6. Core locations required: 
3098 bytes 
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COMMENTS AND RESTRICTIONS 


The computation statements are taken from subroutine CLASEM of the DEB TAP 
program . 

REFERENCES 

Documentation and program listing of subroutine CLASEM of the DEBTAP 
program . 
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Subroutine: Compute Keplerian elements (COMPUT) 
PURPOSE 

To compute Keplerian elements from Cartesian elements. 


METHOD 

The Keplerian elements are computed as follows: 

1. Semimajor axis (ELM(l)) 

R 

a = 

2 - (RV 2 /U) 


2. Eccentricity (ELM(2)) 

°=(^ 2 2 ) 1/2 


3. Inclination (ELM(3)) 


i = tan 


(d 2 + d 2 )' 

-1 \ x y / 


1/2 


4. Longitude of the ascending node (ELM(4)) 


U = tan 


i / D 

-1 / x_ 

-D 


5. Argument of perigee (ELM(5)) 

-1 


oo = tan 


R D 
z 


- tan 


iB D - R D 
y X X y 


9 1/2 

- 1 / < 1 - e > f l 
f„-e 2 
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6. Mean anomaly 


*lJ' 


MEAN ANOMALY = tan 


-l/V 


I- L 


where 


R = (r 2 + R 2 + R 2 ) 
\ x y z / 


V = ( v 2 + V 2 + V 2 ) 
\ x y z / 


1/2 


1/2 


D = R V - R V 
x y z z y 


D 


R V - R V 
Z X X z 


D = R V - R V 
z x y y x 


D 


= ( D : 


2 + D 2 + D 2 
x y z 


1/2 


f = 


RV + R V + R V 
xx y y z z 


Vv] 


f = 1 - — 

2 a 

where (R x , R y , R z ) and (V x , V y , V z ) are position and velocity vectors 
in Cartesian coordinates . The quantity U is the gravity constant of the 
Earth. 
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i . hi wwm 


USAGE 


1. Calling sequence: 


CALL COMPUT (RVEC, VVEC, ELM) 

FORTRAN 

Name 

I/O 

Dimension 

RVEC 

I 

3 

VVEC 

I 

3 

ELM 

O 

6 


Description 
Position vector 
Velocity vector 
Keplerian elements 


2. COMMON blocks used: 


CBSTAT, CNVFCT, ECONS 


3. COMMON block parameters used. 


FORTRAN 

Name 

COMMON 

Name 

I/O 

Dimension 

Description 

KTR 

CBSTAT 

I/O 

400 

Subroutine/entry point counter 

DEGPRD 

CNVFCT 

I 

1 

Degrees per radian 

TWOPI 

CNVFCT 

I 

1 

Pi times two 

EMU 

ECONS 

I 

1 

Gravity constant of Earth in 
kmVsec^ 


4. Subroutines and functions used: 

None 

5. This routine is called by: 

Main routine 

6. Core locations required: 

1600 bytes 


COMMENTS AND RESTRICTIONS 

The computation statements are taken from subroutine OSCMN of the DEBTAP program. 


REFERENCES 

Documentation and program listing of subroutine OSCMN of the DEBTAP program. 
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Subroutine: Convert year, month, day, hour, and minute from floating point 
variables to integer variables (DATOUT) 

PURPOSE 

To convert year, month, day, hour, and minute to integer. Seconds remain in 
floating point for output purpose. 

METHOD 

Truncate floating point variables after calling subroutine SECCAL. 

USAGE 

1. Calling sequence: 


CALL DATOUT (REFDMY, SECIN, IDAY, IMONTH, IYEAR, 


FORTRAN 

Name 

I/O 

IHRMIN, SECOUT) 
Dimension 

Description 

REFDMY 

I 

1 

Date in DDMMYY. 0 

SECIN 

I 

1 

Seconds from the dat 

IDAY 

0 

1 

Day 

IMONTH 

O 

1 

Month 

IYEAR 

O 

1 

Year 

IHRMIN 

O 

2 

1 = hour, 2 = minute 

SECOUT 

O 

1 

Second 



2. COMMON blocks used: 

CBSTAT 

3. COMMON block parameters used: 


j L FORTRAN 

COMMON 



] ■ Name 

Name 

I/O. 

Dimension Description 

KTR 

CBSTAT 

I/O 

400 Subroutine/entry point counter 


.. r.'v 

•Vir 
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4. Subroutines and functions used: 

Entry point CALSEC in SECCAL 

5. This routine is called by: 

Main routine 

6. Core locations required: 

6S2 bytes 

COMMEN TS AND RE STRICTIONS 
None 

REF EREN CES 


Documentation of subroutine/ SECCAL of the DEBTAP program 




Subroutine: Compute osculating elements and mean elements (OSCMN) 

PU RPO SE 

To compute osculating elements at ascending or descending node, mean elements 
and other related parameters. 


METHOD 

The Brouwer mean elements required for the averaging are computed from 
Kepler ian elements at each collected point by calling subroutine CONVRT. 
Subroutine COMPOS computes the Kepierian elements for each point used for 
averaging. The averaging applies to the first three of the six Brouwer mean 
elements, namely, semi-major axis, eccentricity, and inclinafT n. The 


averaging operation is given by 

~ = L a./n 
1=1 7 

where n is the number of points used for averaging and a and a refers to either 


semi-major axis, eccentricity, or inclination. 


1. Calling sequence; 

CALL OSCMN (ELE, TIME. RLON, NCNT, SECTOR) 


FORTRAN 

Name 

I/O 

Dimension 

Description 

ELE 

I 

6 

Kepierian elements at descending node 

TIME 

I 

1 

Time of descending node 

RLON 

I 

1 

Longitude at descending node 

ncnt ; ; 

I 

1 

Number of points used for averaging 

SECTOR 

I 

1 

Time interval in seconds between two 
points 


2. COMMON blocks used: 


CBMNVR, CBSTAT. CNVFCT, DUMP, ECONS, ELEM, MORCON, 
RAS 
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3. COMMON block parameters used: 
FORTRAN COMMON 


Name 

Name 

I/O 

Dimension 

Description 

REPTIM 

CBMNVR 

I 

1 

Repeat cycle length in days 

NREVDY 

CBMNVR 

1 

1 

Revolutions per day 

NUMREV 

CBMNVR 

I 

1 

Number of revolutions per repeat 
cycle 

KTR 

CBSTAT 

I/O 

400 

Subroutine/entry point counter 

DEGPRD 

CNVFCT 

1 

1 

Degrees per radian 

TWOPI 

CNVFCT 

I 

1 

Pi times two 

L 

DUMP 

0 

1 

Intermediate output control 
30 = no intermediate output 
6 = intermediate output 

ERAD 

ECONS 

I 

1 

Radius of Earth in km 

EMU 

ECONS 

1 

1 

Gravity constant of Earth in 

kmVsec 2 

X 

ELEM 

I 

100 

Position vector of a point in 
km 

Y 

ELEM 

I 

100 

Z 

ELEM 

I 

100 

XDOT 

ELEM 

I 

100 

Velocity vector of a point in 
km/ sec 

YDOT 

ELEM 

I 

100 

ZDOT 

ELEM 

I , 

V 100 

WE 

MORCON 

I 

1 

Rotational velocity of the Earth in 
rad/ sec 

SPD 

MORCON 

I 

1 

Seconds per day 

RMN 

RAS 

I/O 

6 

Mean elements 

DOT 

RAS 

o 

3 

Secular terms 

PART 

RAS 

0 

3 

Partial derivatives of secular 
terms 

OVERLP 

RAS 

0 

1 

Overlap error 

RPTCYL 

RAS 

0 

1 

Repeat cycle error 

TIMMLT 

RAS 

o 

1 

Mean local time 
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4. Subroutines and functions used: 

COMPOS, CONVRT 

5. This routine is called by; 

Main routine 

6. Core locations required: 

3.758 bytes 

COMMENTS AND RESTR ICTION S 
None 

R EFERENCE S 

Documentation and program listing of subroutine OSCMN of the DEBTAP program 
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^5ip 






C 


GSCMN 



CdNVRT 

OBTAIN OSCU- 
L ATI NO elemei 
AT DESCEND/* 
NODE 



Compute 

■ 

SET SUMM/N6- 

■ 

SET LOOP 

PERIOD $ 

: 

VARIABLES, 

SP 

COUNTER 

MAXIMUM 

• 

e, * , 

W 

XCAJT - 1 

t/me l/m/T 


TO EERO 




-*0 




Figure 3-2. Subroutine OSCMN Flow Diagram 












Subroutine: Pack or unpack a numerical code (PACK) 


PUR POSE 

1. PACK 

To pack an array of a given base into a single code of base 10. 

2. UNPACK (entry point) 

To unpack a single code of base 10 into an array of a given base. 
METHOD 

Packing and unpacking is done according to: 

„ . ai a. 2 .,n-l 

C = A- + A 0 b + A b + + A b 

12 3 n 

where 

C = single code 
A = array 

b - given base I 

n = numbers in an array 

USAGE 

1. Calling sequence for PACK/UNPACK: 

CALL PACK (IARRAY, NARRAY, IBASE, ICODE) 
CALL UNPACK (ICODE, IBASE, NARRAY, IARRAY) 
FORTRAN I/O 


Name 

PACK 

UNPACK 

Dimension 


IARRAY 

I 

0 

Variable 

Array 

NARRAY 

I 

I 

1 

Numbers in an array 

IBASE 

I 

I 

1 

Given base 

ICODE 

0 

I 

1 

Code 

2. COMMON blocks used: 




CBSTAT 
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3. COMMON block parameters used: 


FORTRAN COMMON 

Nam e __Name_ I/O Dimension Desc riptio n 

KTR CBSTAT I/O 400 Subroutine/entry point counter 

4. Subroutines and functions used: 

None 

5. This routine is called by: 

Main routine, YIELD 

6. Core locations required: 

650 bytes 

COMMENTS AND RESTR ICTIONS 

This routine is used to pack and unpack YYMMDD or HHMMSS into components 
and vice versa. Therefore, the given base is 100 and the number in an array 
is 3. 

REFERENCES 


Subroutine: Read data from cards (READER) 

PUR POSE 

To read data from cards and perform end-of-file test. 
METHOD 

1. Title is read in format 10A8. 

2. Mandatory constants are read in namelist INCARD, 

3. Optional constants are read in formats 8G10 and 3G25.16. 

USAGE 

1. Calling sequence: 

CALL READER 

2. COMMON blocks used: 


CBMNVR, CBSTAT, DUMP, ECONS, INDATA. MORCON 
3. COMMON block parameters used: 

FORTRAN COMMON 

Name 

Name 

I/Q 

Dimension 

Description 

Repeat cycle length in days 

REPTIM 

CBMNVR 

I 

1 

NREVDY 

CBMNVR 

I 

1 

Revolutions per day 

NUMREV 

CBMNVR 

I 

1 

Number of revolutions per 
repeat cycle 

KTR 

CBSTAT 

I/O 

400 

Subroutine/ entry point counter 

L 

DUMP 

o 

1 

Intermediate output control 
30 - no intermediate output 
6 = intermediate output 

erad 

ECONS 

o 

1 

Radius of Earth in km 

EMU 

ECONS 

o 

1 

Gravity constant of Earth in 
km 3 /sec^ 

TITLE 

INDATA 

0 

10 

Title record 
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FORTRAN COMMON 


Name 

Name 

I/O 

Dimension 

Description 

DELMIN 

INDATA 

0 

1 

Tentative interval size between 
points in minutes 

NCNT 

INDATA 

O 

1 

Number of points to be used for 
averaging 

IEPYMD 

INDATA 

O 

1 

Optional epoch in YYMMDD and 

IEPHMS 

INDATA 

O 

1 

HHMMSS 

IWRITE 

INDATA 

O 

1 

Intermediate output indicator 
that defines L in common block 


DUMP 

IWRITE = 0, L = 30 
No intermediate output 
IWRITE * 0, L = 6 
Intermediate output 

4. Subroutines and functions used: 

None 

5. This routine is called by: 

Main routine 

6. Core locations required: 

938 bytes 

COMMENTS AND RE ST RICTI ONS 
None 

REFE RENC ES 

None 
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Subroutine: Shift an array (SHIFT) 


PUR POS E 

To shift an array left by 300 positions and fill the rightmost 300 positions with 
another array. 

ME THOD 

Use a DO-LOOP from 1 to 300 to address each element in the array. 

USAGE 

1. Calling sequence: 

CALL SHIFT 

2. COMMON blocks used: 

CBSTAT, $BLANKCOM 

3. COMMON block parameters used: 


FORTRAN 

Name 

COMMON 

Name 

I/O 

Dimension 

Description 

KTR 

CBSTAT 

I/O 

400 

Subroutine/ entry point counter 

TEMP 

$BLANKCOM 

I/O 

900 

Temporary array 

TAPE 

$blankcom 

I/O 

350 

Ephemeris tape elements 


4. This routine is called by: 

Main routine 

5. Subroutines and functions used: 

None 

6. Core locations required: 

316 bytes 

C OMMENTS A ND RESTRI CTIONS 


None 



REFERENC ES 

None 
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Subroutine: Store vector elements (STORE) 


PUR POSE 

To store position and velocity vectors from temporary arrays and convert their 
units from DUL and DUL/DUT to km and km/sec. 


METHOD 

Temporary array elements are stored into x, y, z, x, y, and z by increasing 
index of the array by 1. 

The conversion factors are: 

1 DUL * 10000 km 
1 DUT » 864 sec 
1 DUL/DUT = 1/0.0864 km/sec 


US AGE 

1. Calling sequence: 

CALL STORE (IPTBG . IPTND. TFT. TIME, IND, NSTEP6) 

FORTRAN 


__Name_ 

I/O 

Dimension 

Description 

IPTBG 

I 

1 

First point to be stored 

IPTND 

I 

1 

Last point to be stored 

TPT 

I/O 

1 

Corresponding time 

TIME 

I 

1 

Time interval between points 

IND 

I/O 

1 

Index to point a temporary 
array element 

NSTEP6 

2. COMMON blocks used: 

I 

1 

Increment between points taken, 
times six 


CBSTAT. ELEM, $B LANK COM 


3. COMMON block parameters used: 
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FORTRAN 

COMMON 



Name 

Name 

I/O 

Dimension 


Description 


KTR 

CBSTAT 

I/O 

400 

X 

ELEM 

0 

100 

Y 

ELEM 

0 

100 

Z 

ELEM 

0 

100 

XDOT 

ELEM 

o 

100 

Y'DOT 

ELEM 

o 

100 

ZDOT 

ELEM 

0 

100 

TEMP 

$B LANK COM 

I 

300 

4. This 

routine is called by: 




Main routine 


5. Subroutines and functions used: 


Subroutine/entry point counter 

Position vector of a point in 
km 

Velocity vector of a point in 
km/sec 

Temporary array 


None 


6. Core locations required: 


596 bytes 


COMMENTS AND RESTRICTI ONS 
None 


REFE RENC ES 

None 
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Subroutine: Unpack time codes (YIELD) 


PURPOSE 


To convert YYMMDD and fraction of a day in seconds into year, month, day, 
hour, minutes, and seconds. 

ME THOD 

1, Unpack YYMMDD into year, month, and day. 

2. Convert seconds into hour, minutes, and seconds. 


USAGE 


1. Calling sequence: 


CALL YIELD (1ND, SECTOR, IEPHY, MEPHM, IEPHD, 


IEPHH, IEPHM, DEPHS) 


FORTRAN 

Name 


Dimension 


SECTOR I 

IEPHY O 

MEPHY O 

IEPHD O 

IEPHH O 

IEPHM O 

DEPHS O 

2. COMMON blocks used: 

None 

3. COMMON block parameters used; 

None 

4. Subroutines and functions used; 

PACK/UNPACK 


Des cription 


Date in YYMMDD 

Fraction of a day in seconds 

Year in 1900's 

Month in numbers 

Day in numbers 

Hour 

Minute 

Second 
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5. This routine is called by: 
Main routine 

C. Core locations required: 
596 bytes 


C OMMENT S AND RESTRICTIONS 
None 


REFERENCES 


None 


3. 2 COMMON BLOCK DESCRIPTIONS 


The ERTS-B Stand-Alone Software Package uses nine COMMON blocks whose 
names, along with the names of routines that use these COMMON blocks, are: 

CBMNVR 

ERRLNG. OSCMN, READER 
CBSTAT 

Main Routine, CALSEC/SECCAL. COMPOS, COMPUT, CONVRT, CROSSP, 
DATOUT, DOTPRD, ERRLNG, LAGRIN/lNTRPX, OSCMN, PACK/UNPACK , 
READER, SHIFT, STORE, SUBCAP, UNITV 
CNVFCT 

Main Routine, CALSEC/SECCAL, COMPOS. COMPUT, CONURT, OSCMN, 
SUBCAP 
DUMP 

Main Routine, Block Data, COMPOS, OSCMN, READER 
ECONS 

Main Routine, Block Data, COMPOS, COMPUT, CONVRT. OSCMN, 
READER, SUB CAP 
ELEM 

Main Routine, OSCMN, STORE 
INDATA 

Main Routine, READER 
MORCON 

Main Routine, Block Data, CONVRT, OSCMN, READER 

RAS 

Main Routine, CONVRT, OSCMN 

Since COMMON blocks CBMNVR, CBSTAT, CNVFCT, ECONS, and RAS are taken 
directly from the DEBTAP program, they contain some dummy variables. The 
following variables describe the individual COMMON blocks: 
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1. COMMON/CBMNVR/REAL8( ), INT4(13) 


Variable 

Name 

Dimension 

Location 
(Real *8 ) 

Description 

B/'.b LON 

1 

1 

Base longitude in degrees 

basmlt 

3 

2 

Base mean local time in hours, minutes, 

DUM1 

19 

5 

and seconds 
Dummy variables 

REPTIM 

1 

24 

Repeat cycle length in days 

DUM2 

31 

25 

Dummy variables 

Variable 
Na.m_e_ 

Dimension 

Location 

(Integer *4) Description 

ID LEI 

7 

1 

Dummy variables 

NREVDY 

1 

8 

Revolutions per day 

IDLE? 

4 

9 

Dummy variables 

NUMREV 

1 

13 

Number of revolutions per repeat cycle 

2. COMMON/ CBSTAT/INT4(400) 


Variable 

Name 

Dimension 

Location 

(Integer *4) Description 

KTR 

400 

' 1 

Sub routine/ entry point counter 


18 = CONVRT 
100 = COMPUT 
151 = OSCMN 
200 = DA TOUT 
203 = CALSEC 
208 = CROSSP 
214 ~ DOTPRD 
233 = LAGRIN 
280 - SUBCAP 
287 = UNITV 
305 = SECCAL 
307 = INTRPX 

350 = READER 

351 “ SHIFT 

352 = STORE 
375 - COMPOS 

398 = PACK 

399 = UNPACK 

400 - ERRLNG 


3. COMMON/ CNVFCT/REAL8 (20) 


Variable 
Name 

Dimension 

Location 
(Real *8 ) 

Description 

DEGPRD 

1 

1 

Degrees per radian 

DUM1 

15 

2 

Dummy variables 

PI 

1 

17 

Pi 

TWOPI 

1 

18 

Pi times two 

DUM2 

2 

19 

Dummy variables 


Data Values for COMMON/CNVFCT/ 


DATA PI /3.141592653589793/D 00/ 

, TWOPI /6. 283 1853071 7 95863D 00 / 

, DEG PRD /5. 72957795130823 D+ 01/ 
4. C0MM0N/DUMP/INT4(1) 


Variable 

Name 

Dimension 

Location 

(Integer *4) Description 

L 

1 

1 

Optional intermediate output indicator, 
defined by IWRITE 

1) IWRITE = 0, L = 30 

No intermediate output 

2) IWRITE * 0, L = 6 

Intermediate output 

5. COMMON/ ECONS/REAL8 (683) 


Variable 

Name 

Dimension 

Location 
(Real *8 ) 

Description 

EMU 

1 

1 

3 2 

Gravity constant of Earth in km /sec 

ERAD 

1 

2 

Radius of Earth in km 

DUM1 

2 

3 

Dummy variables 

ESQ 

1 

5 

Square of the eccentricity of the oblate 
Earth 

DUM2 

676 

6 

Dummy variables 

DERAD 

.:./!■ / 

682 

Dynamic radius of Earth in km 

PORAD 

i 

683 

Polar radius of Earth in km 
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Data Values for COMMON/ CNVFCT/ 


DATA EMU /3. 986012 D + 05/ 

, ERAD /6. 378165 D + 03/ 

, DERAD /B. 378155 D + 03/ 

. PORAD /6. 378784 D + 03/ 


, ESQ /6. 693421622966D - 03/ 
6. COMMON/ ELEM/REAL8 (600) 


Variable 


Location 


Name 

Dimension 

(Real *8) 

Description 

X 

Y 

100 

100 

1 

101 

Position vector of a point in 
km 

Z 

100 

201 

XDOT 

100 

301 

Velocity vector of a point in 
km./ sec 

YDOT 

100 

401 

ZDOT 

100 

501 

7. COMMON/INDATA/REAL8(11) s INT4(4) 

Variable 


Location 


Name 

Dimension 

(Real *8 ) 

Description 

TITLE 

10 

1 

Title of ERTS-B in format 10A8 

DELMIN 

1 

11 

Tentative step size in minutes 

V ariable 


Location 


Name 

Dimension 

(Integer *4) Description 

NCNT 

1 

1 

Number of points to be used for averaging 

XEPYMD 

1 

2 

Optional epoch in YYMMDD and 

IEPHMS 

1 

3 

HHMMSS 

IWRITE 

1 

4 

Intermediate output selector 
1) IWRITE = 0, L = 30 




No intermediate output 




2) IWRITE / 0 ( L = 6 




Intermediate output 


8. COMMON/MORCON/REAL8(5) 


Variable 

Name 

Dimension 

Location 
(Real *8) 

Description 

CJ2 

1 

1 

2nd harmonic constant 


Variable 


Location 


Name 

Dimension 

(Real+8|_ 

Description 


CJ4 

1 

2 

4th harmonic constant 

SPD 

1 

3 

Seconds per day 

WE 

1 

4 

Earth's rotational velocity in rad/sec 

EC CM IN 

1 

5 

Minimum eccentricity acceptable for 
subroutine CONVRT 


Data Values for COMMON/MORCON/ 

DATA C J2/+1 , 0823 D-03/ 

, CJ4/-1. 8000 D- 06/ 

WE/7.292115146 D-05/ 

, SPD/8.64 D+04/ 

, ECCMIN/5. 0 D-04/ 

9. COMMON/RAS/REAL8(3l) 


Variable 

Name 

Dimension 

Location 
(Real *8) 

Description 

OSC 

6 

1 

Osculating elements 

RMN 

6 

7 

Mean elements 

DOT 

3 

10 

Secular terms 

PART 

6 

16 

Partial derivatives of secular terms 

DRL 
DRLA 
DR LI 

1 

1 

1 

22 

23 

24 

Latus Rectum and its derivatives in deg/ 
rev, deg/rev-km, and deg/rev-deg 

OMGDT 

DOMGDA 

DOMGDI 

1 

1 

1 

25 

26 
27 

Omega term and its derivatives in deg/ 
rev, deg/rev-km, and deg/rev-deg 

OVERLP 

1 

28 

Overlap in vertical projection by 
percentage 

RPTCYL 

1 

29 

Repeat cycle error in nautical miles 

TIMMLT 

1 

30 

Mean local time 

RNO 

1 

31 

Mean motion in rad/sec 
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VU. OUTPUT 


Output is cm a single page as shown in Figure 8. 
The following variables are nodal parameters: 


X 


Y 

Position vector 

Z 


XDOT 

YDOT 

ZDON 

Velocity vector 

HT 

Height above surface 

LAT 

Latitude in degrees 

LON 

Longitude in degrees 

V 

Magnitude of velocity 

GAM 

Flight path angle 

AZ 

Azimuth 

SMA 

Semi-major axis 

ECC 

Eccentricity 

AINC 

Inclination 

LAN 

Longitude of descending node in degrees 

APF 

Argument of perigee 

SML 

Semi-latus rectum 

RCA 

Point of closest approach 

APR 

Apo center radius 

PER 

Period 

TA 

True anomaly 

EA 

Eccentric anomaly 

MA 

Mean anomaly 

HGA 

Height of closest approach 
(= APR - radius of Earth) 

HCAN 

HC A/1. 852 in nautical miles 
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APHN APH/l.852 in nautical miles 

The following variables are mean elements: 


SMA 

semi-major axis 

ECC 

Averaged eccentricity 

INC 

inclination 

LAN 

Longitude 

APF 

Perigee 

MA 

Mean anomaly 

DLE 

Semi-latus rectum 

DLRDA 

Derivative of semi-latus rectum with respect 
to semi-major axis 

DLRDI 

Derivative of semi-latus rectum with respect 
to inclination 

DOMGDA 

Derivative of omega with respect to semi- 
major axis 

DOMGDI 

Derivative of omega with respect to inclination 

RPTCYL 

Repeat cycle 

APGH 

Apo center radius 

PRGH 

Peri center radius 

ERRLON 

Error in longitude 

ERRMLT 

Error in mean local time 
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